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ABSTRACT 


Second adiabatic theory is used to derive an asymptotic expression for the 
bounce-averaged rate of kinetic energy change (w) experienced by charged par- 
ticles trapped in electric and magnetic fields which conserve the first two in- 
variants M and J of their motion. The restriction that there exist no electric 
field Ej| , parallel to the magnetic field, reduces {w) by one order in the adiabatic 
expansion parameter m/e from its E ( | ^ 0 value. If in addition the magnetic 
field arises solely from, and convects rigidly with, the arbitrary motion of a 
perfectly conducting magnet, we find that (w) is smaller by yet another order 
in m/e. Our results generalize and confirm the conclusions of Hones and Berge- 
son [l965] who consider E (| = 0 and the special case of a uniformly rotating 


conducting magnetized sphere with non-aligned rotation and magnetization axes. 
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CHARGED PARTICLE ENERGIZATION 
BY AN ARBITRARILY MOVING MAGNET 


INTRODUCTION 

Electric fields arising from the uniform rotation of a magnetized conducting 
sphere have been previously examined as a mechanism for the energization of 
charged particles. [Terletzky , 1946; Davis , 1947, 1948; Hones and Bergeson , 
1965.] These studies have attempted to determine whether such electric fields 
(assuming that a star, the sun, or the earth can be modeled as such a sphere) 
can accelerate particles to cosmic ray and auroral energies or explain the high 
energy populations of planetary magnetospheres. 

Cases where the rotation and magnetic axes are parallel [Davis, 1947, 1948] 
and arbitrarily aligned pTerletzky , 1946; Hones and Bergeson , 1965] have both 
been investigated. If the sphere is assumed to be surrounded by a vacuum, 
large accelerating voltages (>10 4 volts for the earth and higher for the sun and 
stars) have been shown to exist along dipole field lines. 

For such a sphere embedded in a plasma medium, however, the ro- 
tation induced electric field component E i( parallel to the exterior dipole lines 
can be cancelled to a large extent by an electrostatic field arising from small 
plasma charge separations along the field lines. Assuming that Ey = 0 exactly, 
Hones and Bergeson [l965] have shown that (w)/ e for trapped particles vanishes 


through second order in e = m/e: in the magnetic field of a dipole rotating with 
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magnetization and rotation vectors arbitrarily aligned, trapped particles lose as 
much energy over one half of their bounce motion as they gain over the other 
half if K || = 0. 

In this paper we first derive from second adiabatic theory [Northrop , 1963] 
an expression for the (w) of a trapped particle. The leading order contribution 
to (w) arises from the presence of a parallel electric field; when infinite plasma 
conductivity parallel to the B-lines is assumed, precluding the existence of E ( | , 
(W) drops one order in e. We next further restrict the problem by assuming 
that the only source of magnetic fields is an arbitrarily shaped conducting magnet 
whose motion conserves the particle's magnetic moment M and second invariant J, 
but is otherwise arbitrary . The pattern of currents in the magnet's frame is 
assumed invariant, so that the magnetic field moves rigidly. In this case energy 
changes (w) occur on a still slower time scale than when E () =0 alone. 

Our principal conclusion is that the slow variation of (w) found by Hones and 
Bergeson for Eh = 0 and a dipole magnetic source is a particular example of 
a theorem applicable to general motions of arbitrarily shaped magnets. 

A GENERAL EXPRESSION FOR (w) 

The general representation of an electric field E(r, t)in terms of the 
electrostatic potential <P( r, t)and the magnetic vector potential A(r, t) is 


= -V0- 


1 d\ 

c dt 


E 


(1) 
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An alternative form of A is 

A(r, t) = a(r, t) V/3(r, t) + Vr( a, /3, t) . (2) 

In Equation (2) a and /3 are functions of space and time such that magnetic field 
lines are the intersections of surfaces a = const, with surfaces of /? = const. 

The magnetic field B is given in magnitude and direction by 

B(r, t) = Va(r, t)xV/3(r, t) . (3) 

The fact that there exists an infinity of pairs a - (3 for a given B is of no concern, 
for the resultant expression for (w) is independent of the choice of a particular 
a - (3 pair. Given an a - /3 pair, Vr in Equation (2) represents this gauge freedom, 
and since (w) is gauge invariant we set r = 0. 

Once functions a and /3 are chosen, it is readily shown [Northrop , 1963] that 
each magnetic field line retains its a, (3 labels as it moves in time through space 
with the field line velocity 

(dp _ da \ B 

n = (at Va - at V P) * • < 4 > 

A point r in space is labeled at time t by the a - fi values of the magnetic field 
line presently passing through it, together with a measure of arc length s(r, t) 
measured along the field line from an arbitrary reference surface. The set of 


independent variables a, /3, s, t can thus be used instead of the set r, t . 
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Since a and /3are constant on field lines, ya and V/3 lie perpendicular to field lines, 
and 


V(0 + 0) 


(<W) 


is thus a potential for the parallel electric field and is constant on a field 
line segment where E (| = 0. 

In adiabatic theory there is associated with each particle the useful quantity 

K - W + e(0 + 0) , (8) 


where W is the particle kinetic energy and e(0 + 0) is evaluated at the particle 
position. For static fields 0=0 and K is the total energy of the particle moving 


in the potential 0. 
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From Equation (8) we determine the instantaneous time rate of change of 
particle kinetic energy to be 


W - K - e(0 + 0) , (9) 

I 

a dot denoting a time variation as seen by the moving particle. 

Let us denote the average of the arbitrary quantity a( t ) over the bounce 
period T of a given particle by (a) , 


<*> 


T” 1 



dt' a(t') 



(10) 


the s integrals extending over a complete guiding center bounce path, which in 
second adiabatic motion very closely follows a magnetic field line. Averaging 
Equation (9) in this way, we obtain for the bounce averaged rate of change of 
particle energy 

(W> = (K> - e [(<£> + (0>] ■ (11) 


Northrop [l963] has shown 



1 aj(a, /3, K, M, t) , 3 v 
eT dt + ®\ e I 


(12) 
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where J is the particle longitudinal invariant defined, using Equation (8), as a 
function of a, (3, K, M( = W^/B), and t as 

J ( a, / 3 , K, M, t ) - m <j) v (| ds 

= (2m) 1/2 <j> [K - e l(a, / 3 , s, t) - MB(a, /3, s, t)] 1/2 ds . (13) 

In Equation (13) £ and 8 are respectively [<t> + 0 ] (r, t) and B(r, t) expressed in 
terms of the variables a, / 3 , s, and t. The (&(e 3 ) in Equation (12) indicates that 
terms neglected are asymptotically of third order in the smallness parameter e 
Since both K and J contain a factor m, 1/e appears on both sides of Equation (12) 
simply for ordering purposes. 

Taking d/dt of Equation (13) (at constant a, /3, K, and M), we obtain 

if = _T [ e 4v ) + m 

The expression for (w)/e, 



<W> /d A M /• . . , 3 n 

V = vr/ + - <^> + ®( e3 ) 


(15) 


J 
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follows by combining Equations (11), (12), and (14). Noting that 


(0(r, t) + <//(r , t>) - (2(a, 0, s, t>) = 


+ a 


dZ 

da 


• dZ 
P 30 


+ v„ 


ii 

(9s 


» 


(16) 


we alternatively write Equation (15) in the form 


<W> 


M /(9B\ / <92 -52 

e W X da + P 30 + v 


H) + 


(17) 


where in terms of the guiding center drift velocity v D and the field line velocity 

v^, a and 0 are (v D - • Va and (v D - • V/3 respectively. 

Requisites for the existence of the invariants M and J are that time variations 
of B take place on the guiding center drift time scale (an e~ 1 time scale) and that 
E be of order e. Under these conditions dB/dt, a, /3, and2 are C>(e) quantities, 
and the leading order contribution to Equation (17) comes from, the E|j (i.e. 31/3 s) 
term. If Ey / 0, we find using Equation (10) 

~ = _ T { 2 ( a - A s u- t) "2 (a, 0, s^, t)| + C>(e 2 ) (18) 


where s u and so are the upper and lower mirror points for a particle with in- 
variants J and M whose guiding center is on the line a, 0. Even a small amount 
of E || can thus produce significant energy changes on the drift time scale. 
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(W>IN SPECIAL CASES 

For E|| = 0, 2 is independent of s as shown in the preceding section and 
Equation (17) reduces to 


& 


M /c?B\ 
e yt/ 



(19) 


From the relationships 



c <9K(a, 0, J, M, t) 
e d/3 


_c_ dj(a , §, K, M, t) _ 
eT d/3 


(20a) 


,-x = 5 dK(q, J, M, t) 

W e da 


eT d/3 


(20b) 


derived by Northrop [l963], there follow the alternative forms of Equation (19): 


$ 

e 

<W> 


M /aB\ c /k, l\ , , 

T\aT/ + -?i * e H ■ 


M /c)8\ C „ /j, , , , x 

[frrj + a(« 3 ) 


(21a) 


(21b) 


1 being the Jacobian operator. All terms in the forms (19-21) are now C>(e 1 2 ) and 

we conclude that the absence of E ( | slows the energy change by one order. 

When the only source of magnetic fields is a moving conducting magnet, 
and the magnetic field moves rigidly with the magnet, dI3/<?t (q, /3, s, t) = 0: 
to an observer moving so as to always remain at point s on the field line q, (3 , 
no magnetic field changes are apparent. Also, since 2 is independent of s, it is 
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the same at any exterior point on a field line as at the point where the line enters 
the magnet. Therefore, if £ is known everywhere on the surface of the magnet, 
(w) can be obtained at all exterior points from any of the forms (19-21). The 
problem is reduced to determining £ on the surface of the magnet. 

Inside the moving conducting magnet we assume that conduction electrons 
are in equilibrium, i.e. 


E + 


v(r, t ) x B 


1 <9A 

c dt 


v x (y x A) 
c 


( 22 ) 


where v(r, t) is the velocity of a point in the magnet located at r in the inertial 
frame ("laboratory frame") in which (w) is to be evaluated. It will be further 
assumed that only the permanent magnetization contributes to A, which then 
moves rigidly with the magnet. The following considerations follow closely 
those of Backus [l956] but allow for arbitrary motion of the conducting magnet. 

Vector manipulation of Equation (22) leads to 

1 d A 1 

V4> + Jt* = ^ [v(v * A) - A • Vv “ A x (v x v) - v • VA] . (23) 

Here A, 0, v, and r are all measured in the laboratory frame; <JA/dt, the rate of 
change of A seen by an observer fixed in the laboratory frame at r, arises from 
both rotation and translation of the magnet: 


«9A 

dt 


Q(t) x A - v • VA . * 


(24) 
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The instantaneous angular velocity Q is defined in terms of the instantaneous 
rotation rate with respect to the laboratory coordinates of any set of axes fixed 
in the conducting body. Further, 


v - v Q (t) + Q(t) x [r - r 0 (t)] , 


(25) 


where v 0 and r 0 are the velocity and displacement of the origin of the body axes 
with respect to the origin of laboratory coordinates. Substituting from Equations 
(24) and (25) into Equation (23), we obtain after performing the V~ope rations 


V 


(*-^) = ° 


(26) 


This result is identical with that of Backus , who considered the special case 
where the axis of rotation is fixed in the laboratory frame. 

Thus, interior to the conducting magnet 


0 = 


v • A 


+ D(t) , 


(27) 


where D(t) is arbitrary, so that in this region 


<f> + yff = 


v • A 1 d/3 

+ ~T~ + ^ a JT + D(t) 


a 

c 


d/3 

dt + v • V/3 + D(t) 


■] 


(28) 
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Since the field moves rigidly with the magnet, d/3/dt + v • V/3 s 0 and + \Jj is 
simply an arbitrary function of time. 

Taking for boundary conditions the continuity of <t> and <//, we conclude that 
4> +0 is independent of spatial position outside the magnet, since all exterior 
points are connected to the magnet by magnetic field lines and 4> + 0 is s -independent. 
It then follows that di/da = 0, dZ/dfl - 0, and with d'Q/dt = 0, Equation (19) predicts 

= 0 + Cl (e 3 ) . (29) 

Thus the average rate of energy change drops by two orders in e. 

CONCLUDING REMARKS 

We would emphasize that the principal result of this paper— that the bounce 
averaged rate of energy change of trapped particles is of <S(e 3 ) small in the 
adiabatic parameter— is a theorem, proved for the idealized model in which 
E|j = 0 and B arises solely from and moves rigidly with a conducting magnet. 

If the entire magnetic field surrounding the earth were due to the earth's in- 
ternal currents (and E (( = 0) the wobbling field caused by the earth's rotation 
would not affect trapped particle energies. The existence in reality, however, 
of small finite parallel resistivity and concommitant E )( and the presence of non- 
rigidly convected magnet field sources— e.g. ring, magnetopause, and tail 
currents— push us into the range of validity of Equations (18 and 20) and their 
predictions of much more rapid energy changes. 
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